Pressure dependence of the Boson peak in glasses 



V. L. Gurevich 

Solid State Physics Division, A. F. lojje Institute, 194021 Saint Petersburg, Russia 

D. A. Parshin 

Saint Petersburg State Polytechmcal University, 195251 Saint Petersburg, Russia 

H. R. Schober 

Institut fiir Festkorperforschung, Forschungszentrum Jiilich, D-52425 Jiilich, Germany 

(Dated: February 2, 2008) 

The inelastic scattering intensities of glasses and amorphous materials has a maximum at a low 
frequency, the so called Boson peak. Under applied hydrostatic pressure, P, the Boson peak fre- 
quency, ojb, is shifted upwards. We have shown previously that the Boson peak is created as a result 
of a vibrational instability due to the interaction of harmonic quasi localized vibrations (QLV). 
Applying pressure one exerts forces on the QLV. These shift the low frequency part of the excess 
spectrum to higher frequencies. For low pressures we find a shift of the Boson peak linear in P, 
whereas for high pressures the shift is oc P^^^. Our analytics is supported by simulation. The results 
are in agreement with the existing experiments. 

PACS numbers: 61.43.Fs, 63.50+x, 78.30.Ly 



I. INTRODUCTION 

One of the most characteristic properties of glasses is 
a maximum in the low frequency part of their inelastic 
scattering intensities.^ This maximum, the Boson peak 
(BP), originates from a maximum of the ratio g{uj)/uj'^ 
where g{uj) is the density of vibrational states which itself 
often has no corresponding maximum." The BP shows an 
excess of low frequency vibrations above the Debye con- 
tribution of the sound waves. The BP is observed in 
experiments on Raman scattering of light and inelastic 
neutron scattering within the frequency interval 0.5 — 
2 THz. It is considered to be one of the universal prop- 
erties of glasses and is found also in a number of other 
disordered systems (see Ref. Q and references therein). 

Disorder affects the vibrational states differently than 
it affects the electronic ones. The main point is that 
the matrix determining the eigenvalues of the harmonic 
vibrations, the squared frequencies, must satisfy the re- 
quirement of mechanical stability (cf. Refn^). In other 
words, all eigenvalues must be positive, apart from the 
six zero values for translation and rotation of the system 
as a whole. Other than in the case of electronic states, 
there is a fixed zero level. An arbitrary random matrix 
has no such property. This means that, in general, such 
a matrix corresponds to an unstable vibrational system. 
Such a mechanical instability has often been observed in 
numerical simulations of disordered vibrational systems 
for sufficiently large degree of disorder {cf. Refs. [alaQl)- 

However, stability is restored automatically when the 
effects of anharmonicity are taken into account.'' A ran- 
dom matrix with the desired stability property is gener- 
ated in a natural way by solving the corresponding non- 
linear problem. It is remarkable that the density of states 
g{oj) of such a "stable matrix" possesses the BP feature. 
The BP is a reminder of the former mechanical instability 



in the system. 

For the proper interpretation of the BP, the key prob- 
lem is the nature of the vibrations that contribute to 
g{uj). Since the term BP is used for any peak in the low 
frequency inelastic scattering intensity one has to dis- 
tinguish between different cases. In some materials the 
BP is ascribed to low lying optical or transverse acoustic 
modes of parental crystala^iSiifli or to librations of some 
molecules in plastic crystals If these excitations have 
a small frequency spread they will show as a peak in g{Lu). 
The role of disorder is merely to broaden modes which 
exist already without disorder. In the opposite case the 
excess of low frequency modes is caused by disorder itself. 
A simple example is realised in a metallic like model glass 
where the the parent crystal is fcc.^'^ In the present paper 
we discuss this latter case. Between these two extreme 
cases of well defined low frequency modes, broadened by 
disorder, and no such modes before disorder, there is a 
range of materials having aspects of both. Despite being 
derived for the case of a disorder-induced BP, our results 
will apply, at least semi-quantitatively, to these interme- 
diate cases as long as the frequency spread of the relevant 
low frequency modes is sufficiently large. 

Starting point of our investigation are the ubiquitously 
occurring quasilocal (or resonant) harmonic vibrations 
(QLV). These vibrations can be understood as a low fre- 
quency vibration of a small group of atoms which has a 
weak bilinear interaction with the continuum of acoustic 
vibrations of the whole system. They share many prop- 
erties of the localized vibrations but are different from 
these exact harmonic eigenstates. They can be seen as 
resonances in the low frequency part of the local spectra 
of the set of atoms involved. In Mofibauer experiments 
one observes anomalously large Debye- Waller factors for 
the atoms vibrating with such QLV—. In simulations 
the QLV are seen in the harmonic eigenstates as "low 
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frequency localized modes" or mixed into the extended 
modes. This is found for the textbook case of a heavy 
mass defect as well as for QLV in glasse^i^. A "low fre- 
quency localized mode" of a small system does not vanish 
when the system becomes larger - it just hybridizes with 
the other modes of similar frequency. If one looks at the 
modes around the QLV-frequency one still finds an ad- 
ditional eigenmode. A more detailed discussion of QLV 
with emphasis on glasses can be found in our previous 
work on the Boson peak^ where one can also find sup- 
porting references. 

In the soft potential modeliSiii one describes the 
low frequency vibrations of the defect system in the 
harmonic approximation utilizing a basis of (extended) 
sound waves and local oscillators. The latter are the 
cores of the QLV or "bare QLV" which have been found 
to extend over several atomsiSiiS. In this basis there is a 
bilinear interaction which is treated as a perturbation. 



II. EFFECTS OF FINITE CONCENTRATION 
OF QUASI LOCALIZED VIBRATIONS 

In a glass one has a finite concentration of QLV's. As 
discussed in our previous paper"^, this has a profound 
effect on their density of states (DOS). The interaction 
of the QLV's with the sound waves induces an elastic 
dipole interaction between them. First, the interaction 
of soft QLV's with surrounding QLV's of higher frequency 
may lead in harmonic approximation to unstable modes. 
Stability is restored by the anharmonic terms. This leads 
to a linear density of states, g{uj) oc w, for uj < Uc where 
LOc is determined by the typical interaction strength. 

Secondly, these renormalized low frequency QLV's in- 
teract with each other. In effect this means that the 
QLV's are subject to random forces, /, later referred to 
as internal forces. Due to the high susceptibility of low 
frequency vibrations their low frequency DOS is changed 
to g{uj) cx uj^ for oj < ojh- This is a general property of 
the low frequency DOS of non-Goldstone bosonic exci- 
tations in random media™ and has been discussed in nu- 
merous papers on l ow energy excitations in glasses, see 
Refs. 4,17 18 20 21]. It is associated with the so-called 
sea-gull singularity in the distribution of the stiffness con- 
stants of the QLV^^'^® in the soft potential modelii. This 
sea-gull singularity has also been observed in computer 
simulations^^-^''. 

By the crossover between the two limiting regimes of 
the DOS the BP is formed and obtains a "universal" 
shape. The frequency of the BP, Wb, is again determined 
by the interaction strength and thus by the characteristic 
value of the internal forces, /o (in Ref. this quantity 
was denoted as 6f): 

oc /q/^. (1) 

With higher interaction uj^ is shifted upwards and the 
intensity of the BP is reduced^. 



One of the most interesting properties of the Boson 
peak in glasses is its shift towards higher frequencies un- 
der application of hydrostatic pressure. In the present 
paper we will show that such a shift can be visualized 
as a result of a simple physical mechanism. If one ap- 
plies a pressure P onto a specimen of glass the internal 
random forces acquire additional random contributions 
A/ proportional to the pressure. As a result, one ob- 
tains a "blue" shift of the Boson peak frequency with 
pressure. For small pressure this pressure contribution 
is small compared to the characteristic value of internal 
forces /o and the shift of the BP is linear in the pressure. 
With increasing pressure A/ can become larger than /o 
and the BP shifts as Wb oc P^/'^. 

III. RANDOM FORCE DISTRIBUTION IN A 
GLASS UNDER PRESSURE 

In the present Section we will briefly derive, in anal- 
ogy to Ref. 0, the random force distribution under an 
applied hydrostatic pressure. From Hook's law one gets 

e^k = -{P/iK)5,k. (2) 

where Eik is the strain tensor and \/K is the compress- 
ibility of the glass. 

The interaction of a QLV with the strain is bilinear- 
Hint = MkSikX = -{P/3K)Aux (3) 

where Aik is the deformation potential tensor and x is 
the coordinate of the QLV. For simplicity wc will write 
in the following A instead of A^^. Thus the additional 
contribution to the random force due to applied pressure 
is proportional to the pressure 

A/ - {P/3K)A. (4) 

The deformation potential A of a QLV is a random quan- 
tity. In particular, it has a random sign, so that the corre- 
sponding distribution function -D(A) is an even function 
of A. As a result the distribution of the random forces 
/ in the glass remains an even function of / when the 
pressure is applied. 

The total random force / is a sum of two contributions 

/ = / + A/. (5) 

Here / is the internal random force in the absence of 
pressure. If the distribution of the internal forces / is 
Q{f) then the distribution of the total random force / in 
a glass under pressure is given by the convolution 

oo 

Fpif)^ J dAQ(j~^P^D{A). (6) 

— oo 

For P = it reduces to the unperturbed distribution 
Q{f) since the distribution -D(A) is normalized to unity. 
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In the Appendix we give the results of the convolution 
for three cases which should be typical, namely two types 
of Lorentzian distributions and a Gaussian distribution 
of / and A. 



IV. THE BOSON PEAK SHIFT UNDER 
PRESSURE 



Let us discuss how random forces change the frequen- 
cies of the QLV's. In a purely harmonic case, the lin- 
ear forces would not affect the frequencies. Anharmonic- 
ity, however, renormalizes the relevant part of the spec- 
Although the QLV's are essentially harmonic 



tru: 

vibrations their frequencies under applied forces can be 
shifted as in the usual quasi-harmonic approximation. 

To start the description of a QLV, consider an an- 
harmqnic oscillator under the action of a random static 
force / 



Uix) = + MLu(xy2 ~ fx 



(7) 



Here ui is the oscillator frequency in the harmonic ap- 
proximation and A is the constant of anharmonicity. The 
role of the omitted third order term as well as that of the 
distribution of A will be discussed further on. The force 
/ shifts the equilibrium position from a; = to 7^ 0, 
given by 



Axl + MujIxq - / = 0, 



(8) 



where the oscillator has a new harmonic frequency 

c^Lw =^l + 3Axl/M. (9) 

If gi(tJi) is the distribution function of frequencies wi 
and Fp{f) is the distribution of random forces in a glass 
under pressure, then the pressure dependent renormal- 
ized DOS is given by 

00 00 

gp{uj)^ gi{uJi)duJi / djFp{f)6 {u) - ujncw) ■ (10) 



Integrating the (5-function with regard of Eqs. © 
and ^ we get 

= 2mM-^ J~^^^^FpiU^.J. (11) 



J Ju^-uf 
^ 



Here 




(12) 



As shown in Ref. Q and shortly discussed in section 
II the interaction of high and low frequency QLV lead 
to a linear DOS in the relevant frequency range. In 
writing Eq. [7|we took QLV already including this effect. 



Therefore, we can approximate giluji) 
Eq. in the form 



CuJi and write 



<?p(w) = 2CM\l — u; 




xFp 



(1 + 22/2)71^ 



(13) 



Let us introduce a characteristic value fo{P) of the to- 
tal random forces acting on the QLV's under pressure. 
This is the characteristic scale of variation of the argu- 
ment of Fp{f). Depending on the relative strengths of 
the internal forces without pressure and the pressure in- 
duced forces, the following estimates hold 



MP) 



/o [1 + 0(P/Po)] for P«Po 
(P/Po)/o for P»Po. 



(14) 



Here Pq is the characteristic pressure, which in the simple 
case of Lorentzian distributions centered around (see 
Eqs. ED Ell), is given by Pq = S/C/o/Aq. Thus for P < 

Po, /o(P) = fo while for P » Pq, /o(P) cx P. 

According to Ref. for a Lorentzian distribution of 
random forces, the Boson peak frequency is given by 



Wb(P) 



1.9Ai/6/o^/-^(P) 
M172 



(15) 



For frequencies uj <C i^h{P) the argument of the func- 
tion Fp in Eq. H13|) is much smaller than the typical 
value of /o(P). One can replace this function by Pp(0) w 
l//o(P). According to Eq. gp{uj) cx Pp(0)u;4. Thus 
at high pressures P ^ Pq, 



gp{uj) oc iv'^ / P for u uJh{P) . 



(16) 



In the opposite case uj ^ CLib(P) the integral over y is 
dominated by such values of y near the upper limit that 




< 1. 



(17) 



After integration, making use of the normalization of 
Fpif), we regain the equation for the unperturbed 
pressure-independent linear density of states 

gp{Lu) —gi{uj) (X Lu for lu]^(P) <ti uj < tOc (18) 

as it should be. For higher uj {uj > Uc) the linear DOS 
produced by the interaction between the QLV will be 
modified and gp{uj) shows material dependent devia- 
tions. 

Fig. n shows the frequency dependence of the BP for 
different applied pressures. Both the distributions of in- 
ternal forces / and of the deformation potentials A were 
approximated by Lorentzians centered at 0. After convo- 
lution a broadened Lorentzian is obtained, Eq. ljA3|l . One 
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FIG. 1: Boson peak according to Eq. (I13II for different applied 
pressures, excluding the low frequency Debye contribution. 
Lorentzian distributions were assumed for both internal and 
pressure induced forces, Eq. (A3). 



the distributions. The crossover, on the other hand, does 
somewhat depend on the type of distribution used. This 
indicates a material dependence in this pressure range. 
The shift of cjb with pressure can be described for the 
Lorentzian distributions, Eq. (A1-A4), by 

^b(p)=c.b(o)(i+i^y^' (19) 

with 

Po = 3Kf o/Ao. (20) 

and for the Gaussian distributions, Eq. (A5-A8), by the 
slightly different form 

J (21) 

which describes a sharper crossover. 



clearly observes the pronounced flattening with pressure 
of the low frequency part of the BP. Contrary to this, 
the high frequency part is not affected. In calculating 
the curves in Fig. ^ we assumed gi(uj) o: uj. In real mate- 
rials, this linearity of the DOS will only hold up to some 
frequency uj(.- Above that frequency the DOS will be 
material dependent and, therefore, the BP will no longer 
have a universal form above that frequency. This can be 
seen by comparison of the present Fig. ^ with Fig. 3 of 
Ref. 0. 





Eq. (A3) 






- Eq.(A7) 






Eq. (A10) / 












// 
















0.01 


0.1 1 10 

p/p„ 


100 



FIG. 2: Frequency of the Boson peak maximum as a func- 
tion of apphed pressure for different force distributions. Solid 
line: Lorentzian distributions Eq. (A3); dashed line: Gaus- 
sian distributions Eq. (A7); dotted line: double Lorentzian 
distribution, Eq. (AlO), with A = ^Aq. 

The variation of the BP frequency, Wb, with pressure 
is shown in Fig. |21 for different distributions of the in- 
ternal forces and deformation potentials. The limiting 
behavior for small and large pressures is independent of 



V. NUMERICAL SIMULATION 
A. Pressure dependence 

To test our analytic description we extended the nu- 
merical simulations of Ref. 3] to include additional exter- 
nal forces A/, Eq. (HJ. We placed N = 2197 oscillators 
with frequencies < < 1 on a simple cubic lattice with 
lattice constant a — 1 and periodic boundary conditions. 
The bilinear interaction between two oscillators, is 
written as 

f^h?t = 5y(^/4>.x,, (22) 

where Vij is the distance between the oscillators and J 
is the strength of their coupling which results from the 
coupling of bare QLV's and sound waves, J = A? / pv^ . 
Here p is the density and v the average sound velocity of 
the material. To simulate random orientations of the os- 
cillators we took for gij random numbers in the interval 
[—0.5,0.5]. The masses, M^, and anharmonicity parame- 
ters, Ai, were put to 1. The DOS for the noninteracting 
oscillators was taken as goiuj) oc cj", with n = 1, 2. Ran- 
dom forces 

A/^ = g. ■ fL (23) 

were exerted on the oscillators where gi were random 
numbers in the interval [—0.5,0.5]. 

Generalizing the potential energy Eq. |7J for one oscil- 
lator to the system of N oscillators and adding the inter- 
action terms described by Eq. H22I) . we then minimized 
the potential energy of the total system of N coupled 
anharmonic oscillators. In the usual harmonic expansion 
around this minimum we calculated the DOS for different 
values of /extj representing different external pressures. 
This was repeated for up to 10000 representations. 
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The frequency dependence of g{ijj)/uo'^ is given in Fig. 13 
for Qoiijj) oc Lo and J — 0.07, for different strengths of 
the external force, /^^t- The behavior of the analytic 
results, Fig.^ is reproduced. The slightly different max- 
imal intensities for the higher forces (pressures) result 
from the different choice of distribution for the external 
forces (square instead of Lorentzian) . The internal forces 
originate in the simulation directly from the bilinear in- 
teraction. 




FIG. 3: Simulated g{u])/uP for different external force 
strengths (A'' = 2197, go(t^) oc tj, J = 0.07). Curves from 



left to right: /ex, 
10"^ 



10"-*, 6-10' 



10' 



3-10-^ e-lO-" and 



The pressure dependence of BP frequency is shown in a 
double logarithmic plot in Fig.^for two different original 
DOS's: go{i^) lo and go{ijj) oc uP'. The results for 
both sets agree with our theoretical predictions. This 
illustrates that our results are indeed independent of the 
choice of the initial go{uj), as long as it is not too strongly 
peaked. As stressed in our previous work, above some 
frequency lJc which is generally well separated above Wb, 
the redistribution of frequencies becomes ineffective and 
the original go{Lo) survives. For instance in a plot of 
g{Lu)/uj'^, corresponding to Fig. O for .9o('^) oc uP' the 
curves converge to a constant given by the normalization, 
and the maximal intensities decay more slowly. 



B. Oscillator participation numbers 

In our model the formation of the BP is driven by the 
interaction between soft oscillators (bare QLV's). At low 
frequencies this interaction is weak and, therefore, the 
QLV's will be only weakly coupled. The BP frequency is 
determined by the typical interaction strength between 
the oscillators. To quantify the interaction we introduce 
an oscillator participation number 



:H = ( El 



J|4 



(24) 
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FIG. 4: Boson peak frequency, li^b, versus external force 
strength on a double logarithmic scale: crosses: = 2197, 
g,i[uj) Qnup-, J = 0.1; circles: iV = 2197, po(t^) oc lj, J = 0.07. 
The dashed lines give fits with Gaussian force distributions 
(fit parameters u)\, and Pq)- 



where denotes the component on oscillator j of an 

eigenvector of the coupled oscillator system and ( )^ 

indicates the average over all eigenmodes of frequency lo. 
Note that this oscillator participation number is differ- 
ent from the usual (atomic or molecular) participation 
number of an eigenmode of an atomic system. First a 
QLV (an oscillator in the present description) has typi- 
cally an atomic participation number of ten or moreiS*iS. 
An oscillator participation number of ten is then equiv- 
alent to an atomic participation number of a hundred 
or more. Secondly the participation numbers are further 
increased by the interaction between the QLV's and the 
sound waves^S-. Here this hybridisation is only included 
in so far as it brings about the interaction between QLV's. 




FIG. 5: Average oscillator participation numbers as function 
of frequency for different external force strength. The curves 
correspond to the systems of Fig. |3| 



For all values of the applied external force, 



Gxt; "-osc 
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shows in Fig. [S] the same quahtative behavior as func- 
tion of Lo. For small frequencies, one has more or less 
isolated QLV's (nosc(w) ~ 2). With increasing frequen- 
cies, coupling and hence nosc(^) rapidly increases. It 
reaches its maximum around and drops to a plateau 
with nosc(w) ~ 8. This might at first sight look surpris- 
ing, since the coupling between the oscillators, Eq. H22|l 
was not changed. On closer inspection of the coupled 
equations of motion and the equilibrium condition, one 
sees, however, that the external force does in fact change 
the coupling between the single oscillators. The max- 
imal value of nosc(w) increases with opposite to 
3(wb)/w^. This is what one would intuitively expect from 
an increased coupHng. As in the case of the DOS, also 
nosc(i^) depends on the original DOS, go(y>), for frequen- 
cies w > Wc > Wb- 



C. Distribution of anharmonicity parameters 

So far, we have taken the anharmonicity parameter 
A in Eq. 13 as a constant and have neglected possible 
third order terms. To check the influence of distribu- 
tions of these terms, we did additional simulations where 
we introduced distributions of these parameters. The 
scaled results are summarized in Fig. |3 where we also 
show for comparison the theoretical result of Ea.llllfdot- 
ted line). The solid line gives the simulation results for 
go{(-o) oc and J = 0.2 with a fixed value A — \. The 
simulated BP is slightly wider than the theoretical pre- 
diction. The upturn at Lo/ujh ~ 2 indicates the upper 
limit LOc where the interaction strength no longer suffices 
to destroy the assumed density of the non-interacting 
oscillators, go{uj) cx uj"^. Taking, for the same parame- 
ters, random Ai from the interval [0.7,1.3], no significant 
change can be discerned. This is in agreement with our 
previous resulli^ that this anharmonic term provides the 
mechanism to stabilize the interacting oscillators but, it's 
magnitude does not determine the resulting spectrum. 

The situation is different when we add a third order 
term to the energy of the single oscillators, Eq. |7| 



Ui{x) = Ax^/A + BiX^/i + MLufx^/2. 



(25) 



The dashed line in Fig.|Hlshows the resulting BP for Bi — 
hiUJi with bi a random number from [-1,-1-1] (5o('j-') ^ 
Compared to the curves without this term, the BP is 
considerably broadened even though it retains its general 
shape. 



VI. COMPARISON WITH EXPERIMENT 

Unfortunately not too many experimental data are 
available. Our theory should, therefore, be considered 
rather as a prediction concerning future experiments 
than as an interpretation of the existing experimental 
data. A general increase of Wb has been observed in 
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FIG. 6: Boson peak frequency, uj\j, versus external force 
strength on a double logarithmic scale: crosses: N = 2197, 
go{Lu) (XLo'^, J = 0.1; circles: iV = 2197, go{uj) oc u, J = 0.07. 
The dashed lines give the fit with a Gaussian force distribu- 
tion. 



experiments on a number of materials, e.g. 8102^^*^, 
GeOa^^, GeS2^, polybutadiene^^, polystyrene^^ and 
teflon^. Similar shifts have been reported from com- 
puter simulations of Si02^i*^. However most of these 
data are not sufficient for a quantitative analysis. 

The shift of the BP over a large pressure range has 
been measured in a-Si022?. As shown in Fig. [7| the 
experimental data can be fitted by our theory using 
Eqs. H19|) and H20|l assuming a Lorentzian distribution 
[see Eq. (jA3|l ]. The agreement between the theory and 
experiment remains good even for high pressures. Re- 
garding very high pressures, our theory is applicable as 
long as the short range topology that determines the 
structure of QLV's does not change. 
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FIG. 7: The Boson peak in a-Si02 under pressure; filled 
squares are data of Ref. 

Fig. |S1 shows the shift of the Boson peak in a-GeS2 
measured by Raman scattering. Although the scatter of 
the experimental points is rather large, again the general 
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FIG. 8: The Boson peak in a-GeS2 under pressure; open cir- 
cles are data of Ref. ,2^ 



agreement with our theory is encouraging. 
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FIG. 9: The position of the bump Tmax(-P) in the specific heat 
C{T) /T^ in teflon under pressure; open circles are the data 

of Ref. lai 



the Boson peak in glasses and some other disordered sys- 
tems. Anharmonicity stabihzes the structure but does 
not determine the shape of the Boson peak. The vibra- 
tions forming the Boson peak are harmonic. 

The present paper extends these ideas. We show that 
under the action of hydrostatic pressure the Boson peak 
is shifted to higher frequencies. At comparatively low 
pressures the shift is linear in pressure P while for high 
pressures it is proportional to P^^^. These conclusions 
are in good agreement with the existing experimental 
data. Our work explains the shift of the Boson peek 
without the need to postulate additional negative third 
order anharmonicities^. 

To obtain a quantitative proof, more extensive inves- 
tigations of the pressure dependence of the Boson peak 
position in various disordered systems are needed. Since 
the proposed mechanism is very general, it will also be 
interesting to investigate both theoretically and experi- 
mentally the behavior of the Boson peak under different 
types of strain other than the hydrostatic one, studied 
here, as well as under static electric fields. In future 
work, we hope to show that the same physical mech- 
anism is fundamental not only for the formation of the 
Boson peak but also for such seemingly different phenom- 
ena as creation of the two level systems that dominate 
the properties of glasses at low temperatures. 
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Boyer et a/.^° measured the shift with pressure of the 
low temperature maximum C{T)/T'^ in Teflon, where 
C(T) is the specific heat. This maximum is directly re- 
lated to the BBM. a gain the observed shift fits well with 
our predictions, see Fig. El 

However, the experiments on the change of the Bo- 
son peak position under pressure are so far insufficient. 
Therefore, we believe that further detailed investigations 
of this phenomenon are called for. 



VII. CONCLUSION 

In our previous paper- we have proposed a mechanism 
of the Boson peak formation. The essence of the mecha- 
nism is that a vibrational instability of the spectrum of 
weakly interacting QLV is responsible for the origin of 



APPENDIX A: CONVOLUTION 

To make our investigation more general, we will con- 
sider Lorentzian and Gaussian distributions for both the 
random forces and the deformation potential. Besides, 
we will consider also a distribution of the deformation po- 
tential A that may be called "double Lorentzian". This 
is formed by two superimposed Lorentzian distributions 
with widths Ao which are centered at ±A, respectively. 
For A > Ao/a/S the resulting distribution then has two 
symmetric side maxima and a minimum at A = 0. 

For the Lorentzian distributions centered at zero one 
has 

oo 
— oo 
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and 



oo 

= = ^ I rf^2e^^--^°l-l (A2) 



As a result of the convolution of these distributions one 
gets again a Lorentzian centered at zero but with greater 
width 



FpU) 



1 UP) 



where 



/o(P) = /o + 



2>K 



(A3) 



(A4) 



In the same manner the convolution of two Gaussian 
distributions 



QU) 



OO 

2lT J 



(A5) 



and 



D{k) 



27r 



AoV27r 



(A6) 



leads to another Gaussian distribution with increased 
width: 



Fp{f) 



1 



exp[-/V2/2(P)] 



(A7) 



where 



AqP 
3K 



(A8) 



Finally let us convolute the Lorentzian for the internal 
forces, Eq. lAll with a double Lorentzian distribution: 



DiA) 



1 

2^ 



Ao 



Ao 



(A-A)2 + A2 (A + A)2 + A2 
As a result of the convolution one gets 



(A9) 



Fp{f) - 



fojP) 
2tt 



1 



(/ + aA)2 + /2(P) 
1 

(/-«A)2+/2(P) 



(AlO) 



where a = P/3K and /o(P) is given by Eq. (A4). This 
distribution now depends both on the widths of the two 
Lorcntzians and on the distance between their centers, 
2A, in D{A). 
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